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The role of the selection pressure and mutation amplitude on the behavior of a single-species 
population evolving on a two-dimensional lattice, in a periodically changing environment, is studied 
both analytically and numerically. The mean-field level of description allows to highlight the delicate 
interplay between the different time-scale processes in the resulting complex dynamics of the system. 
We clarify the influence of the amplitude and period of the environmental changes on the critical 
value of the selection pressure corresponding to a phase-transition "extinct-alive" of the population. 
However, the intrinsic stochasticity and the dynamically-built in correlations among the individuals, 
as well as the role of the mutation-induced variety in population's evolution are not appropriately 
accounted for. A more refined level of description, which is an individual-based one, has to be 
considered. The inherent fiuctuations do not destroy the phase transition "extinct-alive", and the 
mutation amplitude is strongly influencing the value of the critical selection pressure. The phase 
diagram in the plane of the population's parameters - selection and mutation is discussed as a 
function of the environmental variation characteristics. The differences between a smooth variation 
of the environment and an abrupt, catastrophic change are also addressesd. 

PACS numbers: 87.10.-|-e, 0.2.70.Lq 



I. INTRODUCTION 

The study of population dynamics is a recurrent sub- 
ject of interdisciplinary research ever since the pioneering 
works of Lotka 1] and Volterra Q . Most of the situations 
involve interacting multi- agent systems, that evolve in a 
spatially-extended environment, and exhibit complex co- 
operative behavior, see e.g. [s', Such systems are gen- 
erally far from equilibrium and highly nonlinear. More- 
over, it has been recognized that simple mean-field (MF) 
like approaches are not always able to grasp all the rich- 
ness of the behavior of such systems [5] . The stochastic 
aspect of the dynamics, as well as the use of discrete spa- 
tial and/or temporal variables to characterize the popu- 
lation evolution may play an important role [1, 0|- Also, 
the inherent discrete nature of a population, i.e., the fact 
that it consists of a finite, integer number of individuals, 
was found to be sometimes an influencial factor of the 
dynamics, see [1, |^ . 

Therefore, the methods developed in nonequilibrium 
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statistical physics are particularly well-suited for the 
study of such systems. In particular, the so-called 
individual-based model (IBM) approach is able to account 
for all the above-cited ingredients of the dynamics. 

A major part of the existing literature is investigat- 
ing various evolution mechanisms in a constant environ- 
ment, like, e.g.. Darwinian selection leading to the elim- 
ination of the ill-adapted individuals [To| . genetic her- 
itage and accumulation of diseases mutations p^ . 
development of hunting/avoiding strategies of preda- 
tors/preys |13| , topology of the network connecting the 
individuals [14 1, etc. As a result of the nonlinear stochas- 
tic dynamics associated with the above mechanisms, non- 
trivial collective properties are emerging. In particular, 
depending on the values of the control parameters charac- 
terizing a system, different steady states are possible 
and a phase diagram can be drawn. 

However, a fundamental question concerning the dy- 
namics of populations or ecosystems is their robustness 
with respect to external perturbations, like, for example, 
a change in their environment. Such an environmental 
change can be abrupt (due, e.g., to some catastrophic 
events) or rather smooth and periodic (due, e.g., to slow 
climatic modifications). Although the climatic changes 
recorded in paleoclimatic data could hardly be called 
perfect oscillations, one notices a certain degree of peri- 
odicity. Periodic changes of the habitat are also assumed 
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in most of the studies done both by biologists and physi- 
cists, see e.g. [l3,[H,[l3], and can be regarded as a rather 
satisfactory approximation to the real data, as far as ex- 
tracting essential qualitative features is concerned. 

One expects on general grounds that new complex be- 
havior of the populations will emerge due to the com- 
petition between the characteristic intrinsic time scales 
and nonlinearities of the system and the time scale of 
the environmental perturbation. This is an aspect that, 
to our knowledge, has not been investigated systemati- 
cally in the existing literature. Indeed, most of the pre- 
vious studies on related points (see, e.g., [1, H^jHH) are 
only reporting on the observed complex dynamics, with- 
out further analysis of its rooting. In this paper we shall 
thus address the role of selection pressure and mutation 
amplitude on the relevant properties of a single-species 
population which lives in a periodically changing habi- 
tat. Their interplay with the magnitude and period of 
the environmental changes will be a fundamental point 
of the analysis. 

The paper is organized as follows: in the following sec- 
tion we present the details of the model; its implemen- 
tation, both in a MF-like and an IBM-type of approach, 
are described in Sec. IIIII Sec. IV presents the analytical 
and numerical results obtained in the MF approxima- 
tion. Section V resumes the main results of the IBM 
simulations, and compares them to the results of the MF 
approach. Finally, conclusions and perspectives are rele- 
gated to Sec. VI. 



II. THE MODEL 

We consider a population that consists of 
hermaphrodite individuals living in a spatially- extended 
habitat, which is represented here as a two-dimensional 
square lattice of size L x L. We assume non-permeable 
boundary conditions, meaning that the individuals 
cannot cross the borders of the lattice. Moreover, the 
lattice has a finite carrying capacity, which comes from 
an exclusion assumption, according to which there is at 
most one individual in each lattice node. 

The dynamics of the population is the result of several 
elements, namely natural selection, individual motion, 
mating and reproduction, whose modeling is described 
below. 

A. Natural selection. Individual trait, optimum, fitness, 
selection pressure, extinction probability. 
Each individual is characterized by its trait, or pheno- 
type, which may correspond to various specific proper- 
ties, such as e.g. the metabolic rate, the body size, or 
may encode a rather complex ensemble of properties. 
The trait of the "i-th" individual is represented here, 
for simplicity, through a real number Zi G [0,1]. This 
number is fixed at the individual's birth, and remains 
constant. This simplified way to characterize an individ- 
ual has been often used both by biologists [l^, [2l| and 



by physicists [25|. Its comparison with more elaborated 
descriptions, that take into account further details of the 
genetic structure of the individuals, was addressed, for 
example, in [24j. 

The population lives in an environment whose influ- 
ence on the individuals is encoded in the value of another 
number ip S [0, 1], the so-called optimum. The optimum 
is that value of the trait which guarantees for its owner 
maximum chance of survival. In principle, the optimum 
could be spatial and/or time-dependent. In this work the 
optimum has the same value in all the sites of the lattice, 
however its value can vary in time. 

The degree of agreement of the individual trait Zi with 
the optimum determines the individual extinction proba- 
bility per unit time (or extinction rate), 
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where 5 is a parameter which models the selection pres- 
sure (SP) of the environment and constitutes a main con- 
trol parameter of the system. The fitness fi of the indi- 
vidual i is defined as 



fi = i- \zi - (p\ ■ 



(2) 



A perfectly- adapted individual has a fitness equal to 1, 
and thus a minimum possible extinction rate. A fitness 
smaller than 1 corresponds to a worse adaptation; fi goes 
to zero for the completely unadapted individuals, that 
have the maximum possible extinction rate, equal to po- 
A variation in time of the fitness of one individual is 
resulting only from the variation of the optimum. 

The constant po in the definition ([1]) of the extinc- 
tion rate is related to the choice of the unit time and it 
depends on whether we are considering the continuous- 
time MF-type model, or the discrete-time IBM-type of 
approach, see below. 

The choice we made of the extinction rate ([!]), the 
implicit definition of the SP parameter S and that of the 
fitness ([2]) are frequently encountered in the biological 
literature, see e.g. [2l|. Other choices, and thus other 
ways of measuring the "selection pressure" and the 
"individual fitness" , are possible. However, most of 
them can be mapped one onto the other and/or account 
for equivalent qualitative aspects of the interaction 
between the individuals and their environment. 

B. Individual motion. 

If an individual survives, it can move to its surroundings. 
The simplest possibility, that we shall adopt hereafter, 
is a random-walk, diffusive-like motion. For example, 
in the discrete-time IBM-type of approach, in one time 
step the individual jumps on the lattice, from its initial 
location to a randomly chosen nearest-neighbor one (i.e., 
a site within the von Neumann neighborhood of the 
initial node) , provided that the chosen site is empty, and 
that it lies within the boundaries of the system. If none 
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of the four first-neighbor nodes is empty, then the in- 
dividual cannot move, and thus cannot mate (see below). 

C. Mating and reproduction. Heredity and mutation. 
Suppose an individual i reaches a destination node. 
Then, if there are other individuals ( "neighbors" ) in the 
nearest-neighborhood of this destination site, the indi- 
vidual "i" choses at random one of these neighbors (call 
it "j") for mating [^. The pair of individuals i and 
j may afterwards give birth to offsprings, whose number 
cannot exceed a prescribed value iVoff . At their birth, the 
progenies are placed, at random, on the empty nodes of 
the joint nearest-neighborhoods of the two parents (that 
counts 6 sites); therefore, the maximum number of off- 
springs A'off ^ 6. If there is no room in this neighborhood 
for putting an offspring, then this one is not born. 

The trait of a progeny k coming from parents i and j is 
determined by the parents' traits (heredity), but it can 
also present some "variations" due to different random 
factors, such as recombination, mutations, etc. We shall 
consider 
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where represents these variations. It brings diversi- 
fication into the phenotypic pool of the population and 
we call it conventionally mutation. For simplicity, we 
shall admit that is a random number, uniformly dis- 
tributed in the interval [—M, M], where < < 1 
is called hereafter the mutation amp litude (MA) and is a 
control parameter of the system [3^. Moreover, if Eq. ^ 
leads to Zfe > 1 or < 0, then one renormalizes it by 
resetting it to = 1, 0, respectively. This means simply 
that the trait of the individuals cannot overcome some 
fixed limits. This choice ^ for the trait of an offspring 
is often made in the biological literature [21i] . 

The population dynamics is thus driven by two main 
"forces" that are acting, to some extent, in opposite di- 
rections: selection and mutation, characterized, respec- 
tively, through the values of the control parameters S and 
M. Selection, combined with heredity, tries to bring the 
average trait close to the instantaneous optimum, while 
mutation introduces diversity in the individual traits, 
and thus is broadening the distribution of the popula- 
tion's traits. 



III. IMPLEMENTATIONS OF THE MODEL 

The study of this system can be made on different lev- 
els of modelling. The simplest one, which has the advan- 
tage to allow, up to a certain extent, for analytic investi- 
gations, is a MF-like level. The drawback of this simple 
approach is to neglect fluctuations, which can play a cru- 
cial role. We are therefore proposing an investigation 
based on IBM-type simulations, for which the fluctua- 
tions are naturally built-in. A comparison of the results 
of these two approaches is then made. 



A. Mean-field like approach 

The main quantity of interest in this approach, that 
neglects the possible correlations between lattice sites, is 
the average population number density c{t). One writes 
down a continuous-time evolution equation for c(i) by 
taking into account the different dynamical processes de- 
scribed above, 



dc{t) 
dt 



= c[i?i(c)i?2(c) 
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This represents the balance between births, described 
by the nonlinear term [c Ri{c)R2{c)], and deaths, cor- 
responding to the (— cp) term. The meaning of these 
terms is as following: 

(a) p is the mean extinction probability per unit time, 
and, as resulting from Eqs. ([T|), ([2]), it reads 
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where (z) is the mean value of the trait. The time scale 
will be chosen such that the constant po = 1. 

In agreement with Eq. ([3]), the MF approach does not 
allow for a self-consistent evaluation of the mean trait 
(z); it remains thus an arbitrary constant e [0, 1]. The 
role of the mutation is practically eliminated in the MF 
description, and the only remaining control parameter is 
the selection S. 

(b) i?i(c) expresses the probability that an individ- 
ual jumps, per unit time, to a randomly-chosen empty 
nearest-neighbor site. 



Ri{c) ^ ro {I - c^) 



(6) 



where (1 — c^) is simply the probability that at least one 
of these four neighbor sites is empty. The constant rg 
depends on the choice of the unit time. If the unit time 
is fixed by chosing po = 1, then ro is also fixed, and it 
measures the ratio between the mean survival time of an 
individual and its characteristic diffusion time on the lat- 
tice. 

(c) The term i?2(c) describes the probability that the 
considered individual encounters at least one nearest- 
neighbor at the destination site and it produces a certain 
number of offsprings according to the model rules: 



i?2(c) =iVoffc(l-c) 



(7) 



One is supposing implicitely that the time scale related 
to mating and reproduction is much shorter than both 
the time scales of survival and diffusion of an individual. 

It can be easily shown that rg ^ 1 implies a rapid 
increase of c till the stationary saturation value c = I, 
while ro ^ 1 leads to an extinction of the population. 
In the foregoing we shall thus focus on the case when 
diffusion and extinction have comparable characteristic 
time scales, and we shall fix throughout rg = 1. 
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One is then led to the following MF equation: 

dc{t) 



dt 



7Voff'C^(l-c4)(l-c)^° 
1 — exp 



l-l^-(z) 



(8) 



This one can be mapped onto the equation of motion 
of an overdamped particle in a potential Usci{c) whose 
actual profile depends on the selection parameter 5, 



dc{t) dUsciic) 



dt 



dc 



(9) 



This point of view will prove to be very useful in the 
discussion of the results presented in Sec. IIVI 
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B. IBM- type simulations 

The IBM simulation algorithm considers the individu- 
als distributed on the lattice nodes, the initial condition 
being represented by their positions and the prescribed 
values of the individual traits. Different initial condi- 
tions were considered, see Sec. |V] below. The individuals 
are evolving, at discrete Monte-Carlo time steps (MCS, 
defined hereafter), according to the stages ^-C of the 
dynamics as described in the Introduction, namely: 
^. At a given time t an individual i is picked at random, 
and its extinction probability pi, corresponding to one 
MCS, is determined according to Eqs. ([U, ([H). The con- 
stant po in Eq. ([T]) is chosen equal to 1. Then a random 
number r is extracted from an uniform distribution in 
the range [0, 1]; if r < Pi, the individual dies, otherwise 
it survives. 

B. If it survives, the individual i jumps at random to one 
of the empty nearest-neighbor nodes on the lattice. 

C. Then it possibly mates and produces offsprings. 

If at the time t there are N{t) individuals in the system, 
then the above steps A-C are repeated N{t) times; this 
constitutes one MCS, the unit-time of the simulations. 
Finally, the time is advanced by one step, t ^ t + 1, and 
the above algorithm is repeated. 

In the next two Sections we shall present and then 
compare the results obtained at each of these two levels 
of modeling. 

IV. MEAN-FIELD TYPE ANALYSIS 

A. Constant environment 

Consider first the case of a habitat with constant op- 
timum If = 0.5, and a population described by a mean 
trait (z) — 0.5 [sj. We shall suppose throughout that 
A^ofi = 6. 

One can use simple arguments of the theory of deter- 
ministic dynamical systems to describe the evolution of 
the mean population concentration c{t). The stationary 



FIG. 1: MF stationary solutions for the concentration c as 
functions of the selection parameter S. Continuous lines 
represent the stable solutions c = and c — 0, and the 
dashed line represents the unstable solution c — c'^ . At 
S = Sc ~ 0.417 the stable and unstable solutions dis- 
appear, and for S > Sc the only remaining stationary solu- 
tion is the absorbing state c — 0. The other parameters are 
If — (z) = 0.5, and A^'off = 6. 



concentrations as functions of the unique control param- 
eter S are represented in Fig.[TJ At a critical value of the 
selection S = Sc « 0.417 a discontinuous phase transi- 
tion takes place in the system. For S < Sc, the system 
presents two stable states, namely ^ and the absorb- 
ing state c = 0, and their respective basins of attraction 
are separated by the unstable solution c^. At 5 = 5c 
the stable and unstable critical points of the over- 
damped dynamics (|4]) disappear. For large S > Sc the 
population gets extinct in the long-time limit, no matter 
what was its initial concentration c(0). 

The standard linear stability analysis [26[ allows also 
to determine the characteristic times for the linear re- 
laxation to the stable point , (denoted Tf^^) and, for 
the escape from the immediate vicinity of the unstable 
point , denoted by T^^^. One can also get the analytic 
expression for the linear relaxation time to the absorb- 
ing state c = 0, namely Tj,°[ = [1 — exp(— These 
characteristic times offer an indication on the order of 
magnitude of the nonlinear relaxation times that corre- 
spond to the fully nonlinear dynamics ([9]). Moreover, 
as shown numerically, for a wide domain of intermedi- 
ate values of the selecion pressure < S < Sc, they are 
practically insensitive to changes in S. This key aspect 
is implicitely taken into account in discussing the cases 
(i) and (ii) in Sec. IVB. 

Another point of interest for the non-trivial stable sta- 
tionary solution is its dependence on the maximum 
number of offsprings iVoff . One can retain two important 
aspects, namely: 

(a) At low concentrations (where the exclusion hypoth- 
esis is practically irrelevant), the survival in a highly- 
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demanding environment (large S) is more efSciently en- 
sured when more descendants are born at each mating, 
(b) In a medium- and low-demanding environment, 
is increasing with decreasing Nog (for a given value of 
S). Indeed, a population that gives birth to fewer de- 
scendants at each mating has to compensate the death 
rate (that is fixed by S) through more mating, and thus 
through a larger concentration. However, this concentra- 
tion cannot exceed a certain value, since, on one hand, 
the individuals have to move in order to mate (and the 
corresponding jump probability i?i(c) is decreasing as 
1 ~ c'^), and, on the other hand, there has to be "enough 
room" for the born offsprings to be put-down (indeed, 
the corresponding probability i?2(c) is rapidly-decaying 
as w (1 — c)^°" for large values of c). The stationary 
value results thus from the balance between these con- 
tradictory tendencies. 



B. Time-periodic habitat 

Suppose that the environment is changing periodically 
in time, and consider, as the simplest possibility: 



ip(t) = 0.5-1- A sin 27r 



t - ti 

~t7 



eit - u) 



(10) 



Here A denotes the amplitude of the environmental per- 
turbation (0 < A ^ 0.5), Tp is its period, and is the 
moment of onset of the perturbation; is the Heaviside 
step function. The mean trait of the population is still 
supposed constant, (z) = 0.5, and the selection has a 
certain given value S = Si. 

This situation is formally equivalent to that of a popu- 
lation of constant mean trait in a constant environment, 
{z) — tp ~ 0.5, but for which the S oscillates, with a 
period Tp/2, as: 



S{t) 
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e{t-ti). 
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This point of view allows us to explain easily the behavior 
of the concentration depending on the value of Si and on 
the parameters A and Tp of the perturbation. 

For Si ^ Sc, the population dies after a time of the 
order of T^^^, so this case is not particularly interesting. 

Let us suppose that Si < Sc, and that c(0) is in 
the basin of attraction of the stable solution cf that 
corresponds to Si. After a time of the order of Tf^^, 
the concentration relaxes to cf . We suppose that the 
perturbation of the selection parameter is applied after 
reaching this state, i.e., U ^ Tf^^. After that, S{t) 
starts to grow, and at ti + Tp/A it is reaching the final 
maximum value S{ — Si/ {I — A); then it decreases 
during Tp/A till reaching the initial value tSi, etc. Two 
situations may appear: 



0.27 h 
5 

0.25 



0.23 



0.21 

5 



0.19 



T ^ 100 

— T =1000 

p 



t 1200 1400 1600 1800 2000 

' t 




1100 1150 
t 



1250 



FIG. 2: MF: Possible temporal evolutions of the concentration 
c{t) in a periodically-changing habitat for different values of A 
and of Tp. The other parameters are c(0) = 0.7, Si = 0.3 and 
ti — 1000. The upper panel refers to the case (i) described 
in the text and is obtained for A — 0.2. The lower panel 
corresponds to the case (ii) in the text, and is obtained for 
A = 0.35. 



Case (i): S[ < Sc. 

The amplitude A of the perturbation is such that S{t) 
never overcomes Sc. Then the nontrivial attractor 
(with ^ cf) never ceases to exist, and c always 
remains in its basin of attraction. The population never 
dies, and its concentration is oscillating in time, with 
a period Tp/2 and an amplitude that is generally de- 
termined by both A and Tp, see the upper panel of Fig. [51 

Case (i.a). If the oscillations are very slow, Tp ^ Tf^^ 
for all the values of iSi ^ 5 ^ <Sf , then at each moment 
the concentration relaxes to the instantaneous value of 
C^, see the upper panel of Fig. O Thus, the concentra- 
tion oscillates between the two extremum possible values 
cf ^ c{t) ^ cf (where cff are the equilibrium values 
corresponding to iSij). 

This situation is also illustrated in the second panel of 
Fig. [3] in terms of the potential C/sci(c) that is driving the 
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evolution of c{t), see Eq. ([9]). When S{t) is oscillating 
between Si and iSf, the profile of the potential C/sci(c) 
is also changing periodically. For a very slow evolution 
of C/sci(c), the representative point of the concentration 
c{t) is practically always "stuck" to the minimum of the 
instantaneous potential. 

Case (i.b)li the oscillations of the perturbation are rapid, 
i.e., Tp is comparable to T^.^^, then the concentration does 
not have time to relax to the instantaneous value of the 
attractor , and the amplitude of its oscillation is there- 
fore smaller than in case (a), see the upper panel of Fig. [51 
This amplitude is strongly dependent on Tp, and for very 
small values of Tp the concentration reaches practically a 
constant value. Thus, for a very rapidly-changing envi- 
ronment, one approaches a behavior that is reminiscent 
of the one appearing in the limit of a constant environ- 
ment A = 0. 

The third panel of Fig. [3] represents this situation in 
terms of the competition between, on one hand, the 
relaxation of the representative point c{t) to the instan- 
taneous minimum of [/soi(c), and on the other hand, the 
rapid variation of the position of this minimum: c{t) lies 
always "behind" the minimum. 

Case (ii): Sf > Sc ■ 

Suppose that the amplitude A of the perturbation 
allows the selection to overcome the critical value, 
i.e., Sf ^ S{t) > Sc during a certain interval of time 
T. Then the attractor ceases to exist during this 
interval, and the concentration tends to relax to the 
absorbing state c = 0, on a time scale of the order of T^^^. 

Case (ii.a) If r ^ T^'^j (the case of a slow perturbation) , 
then the population will die during the interval r. 
The corresponding evolution of the concentration is 
represented in the lower panel of Fig. [21 see also the 
upper and middle panels of Fig. [3| for an illustration. 

Case (ii.h). If, however, r < T^^^ (rapid perturbation), 
then the population may survive till the end of the inter- 
val T. At this moment, the attractor starts to exist 
again; if not already out of its basin of attraction, the 
concentration tends to relax to it, till S reaches again 
the initial value iSi and the cycle starts all over again. 
The behavior of the population in this case may be ei- 
ther oscillating (after a transient regime), or may present 
damped oscillations to the absorbing state c = (the lat- 
ter regime, however, is appearing for extremely narrow 
intervals of the parameters). This is illustrated in the 
lower panel of Fig. [2| and also schematically in the upper 
and lower panels of Fig. [3) 

Another important characteristic of the system is the 
extinction time (ET), which represents the time it takes 
for the population to get extinct, i.e., to reach the ab- 
sorbing c = state. The phase transition "extinct-alive" 
is reflected in the behavior of the ET, which is finite for 
SPs below the critical value and becomes infinite at the 
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FIG. 3: MP: The upper panel describes schematically the tra- 
jectory of the system in the plane concentration-selection for 
the case (i) as discussed in the text, for Si = 0.3 and A — 0.2 
(with tfi = (z) = 0.5). The continuous line corresponds to 
the case (i.a) (slow perturbation), while the dashed line rep- 
resents case (i.b) (rapid perturbation). Black arrows refer 
to the increase of <S, and the white ones to the decrease of S. 
The middle panel illustrates the behavior of the representative 
point c(t) (the big black dot) in terms of the oscillating po- 
tential Usisiic) for case (i.a), and the lower panel corresponds 
to case (i.b). The thick arrows sketch the sense of evolution 
of Ubci(c)', the values of 5 corresponding to the represented 
profiles of Usei (c) are also indicated. 
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FIG. 4: MF: Upper panel: schematic representation of the 
trajectory of the system in the plane concentration-selection 
for case (ii) of the main text, for Si = 0.3 and A — 0.3 (with 
if — {z) — 0.5). The continuous line corresponds to case (ii.a) 
(slow perturbation), and the dashed line represents case (ii.b) 
(rapid perturbation). Black/white arrows correspond to the 
increase/decrease of S. Middle and low panels: the behavior 
of the representative point c{t) (big black dot) when the profile 
of the potential Usc\(c) is evolving in the sense sketched by the 
thick black arrows. Middle panel: case (ii.a). Lower panel: 
case (ii.b). The values of S corresponding to the represented 
profiles of f/sci(c) are also indicated. 



critical value of Si, as illustrated in Fig.[5l 

For a constant environment, the ET is finite as long 
as 5i > iSc ~ 0.417, but it tends to infinity for iSj < Sc 
(provided, of course, that the initial concentration c(0) is 
not in the basin of attraction of the absorbing state). 

For the case of a time-periodic environment, the ET 
may also be infinite, or finite. In the first case the con- 
centration c{t) has an oscillatory behavior, while in the 
latter it tends to the absorbing c — Q state. The critical 
value of Si that separates these two situations depends 
in general on both A and Tp, and it is always comprised 
between iSc(l — A) and Sc- For very slow perturbations, 
the critical value of Si approaches iSc(l — A), which is 
the one that corresponds to a constant environment with 
tf = 0.5 -|- A. Note, however, that increasing Tp to the 
limit Tp — > oo (for a fixed A) does not imply that the gen- 
eral behavior of the system "converges" to the one cor- 
responding to a constant environment with Lp = 0.5 + A. 
This point is illustrated in Fig. O the profiles of ET cor- 
responding to increasing Tp values do not approach the 
profile for a constant environment with = 0.5+A = 0.8. 

For rapid fluctuations the critical value of Si depends 
strongly on A and Tp. However, in the Tp — > limit 
(infinitely-fast fluctuations), it approaches the value iSc- 
More generally, the very fast oscillations of the optimum 
(at arbitrary but fixed A) are "smeared out" , simply be- 
cause the relaxational dynamics of the system cannot 
follow these optimum fluctuations. Thus the behavior 
of the system converges to the one corresponding to the 
constant environment if = 0.5. 

For comparison, we also included in Fig. O the ET for 
the case when the optimum presents an abrupt jump from 
tf = 0.5 Sit t < ti to f — 0.5 + A at t > ti, a situation that 
could model a catastrophic change in the environment. 
The interpretation of the corresponding curve "///" is 
obvious. It shows, for example, that a population that 
could survive a rather rapid but smooth variation of the 
optimum cannot survive a catastrophe. 

This detailed discussion clearly illustrates the delicate 
interplay between the time-scale and amplitude of the 
perturbation, on one hand, and the intrinsic time-scales 
and nonlinearities of the unperturbed dynamics, on the 
other hand. This interplay is determinant even in this 
simplified MF-type of approach; we are thus expecting 
it to play a prominent role for the case of a more mi- 
croscopic, individual-based modeling, for which the fiuc- 
tuations (in particular, the mutations) are enriching the 
dynamics even further. 



V. IBM-TYPE SIMULATIONS 

Let us turn now to the results obtained through numer- 
ical simulations in the IBM-type of approach described 
in Sec. IIII Bl One is apriori entitled to expect a richer 
dynamics in this microscopic approach, since: (a) on one 
hand, the mutation is definitely playing a role, and JV[ 
represents the second control parameter besides the se- 
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FIG. 5: MF: Extinction time tcx as a function of the initial 
selection parameter Si for different periods Tp of the pertur- 
bation. The amplitude of the perturbation is fixed to A = 0.3. 
The critical values of Si are comprised between the extremal 
points Sc{l - A) M 0.292 and 5c « 0.417 (represented by 
the big dots on the abcise). The continuous thick curve la- 
belled / gives the ET as a function of Si for the case of a 
constant environment with ip = 0.5, while the dashed-dotted 
curve labelled II is obtained for a constant environment with 
1/7 = 0.8. Finally, the thick dashed curve labelled III corre- 
sponds to an abrupt jump in the optimum, from ip — 0.5 for 
t < ti to LP = 0.5 + A = 0.8 a.t t > ti. The other parameters 
are (z) = 0.5, ti = 1000, and c(0) = 0.7. 



lection S; (b) on the other hand, the fluctuations and the 
finite size of the system can influence in a nontrivial way 
the global dynamics of this highly-nonlinear system. 

The initial conditions for the system, as already men- 
tioned in Sec. IIII Bi are represented by the positions of 
the individuals on the nodes of the lattice, as well as by 
their individual traits. If not explicitly stated otherwise, 
we shall consider hereafter that initially the A^(0) individ- 
uals are randomly, uniformly distributed on the lattice, 
with a given initial concentration c(0) = N{0)/L'^; also, 
their traits are randomly, uniformly distributed in the 
interval [0, 1]. The case when the A^(0) individuals have 
the same trait was also studied. The initial concentration 
is fixed throughout to the rather larg value c(0) — 0.7. 
The role of c(0) on the survival of the population, and the 
existence of a critical initial concentration (the so-called 
minimum viable population concentration) below which 
the population gets extinct for given environmental con- 
ditions are well-known issues, see e.g. [l^HO, [HI, and we 
shall not address them further here. In the foregoing we 
shall keep NoS — 6. 

The control parameters that determine the dynamics 
of the system are the SP - iS and the MA - M.. In the case 
of a variable environment, as described by Eq. pO]) . one 
has to consider two further control parameters, namely 
the amplitude A and period Tp of the perturbation of 
the optimum. The onset time of the optimum's variation 



is fixed to ti = 1000 MCS, and it is sufficiently large 
such that the system relaxes to a stationary state in the 
interval between t = and t = ti. 

Even when the above parameters, as well as the initial 
concentration and mean trait are fixed, there are still var- 
ious possible realizations of the system's evolution, that 
are induced by the different sources of stochasticity which 
are present in the system, namely, (a) the randomness in 
the initial conditions (random positions of the individu- 
als on the lattice, and/or to their randomly-distributed 
individual traits); (b) the stochastic aspect in each in- 
dividual's death; (c) the random, diffusion-like jumps of 
the individuals on the lattice; (d) the mating and the 
stochasticity elements in the progeny-birth, including the 
mutation in the trait of the resulting offsprings. 

Some measurable quantities [(i) to (iv) below] that can 
characterize in an efficient way the behavior of the pop- 
ulations can be cast into two categories: those referring 
to a single realization of the population's evolution, and 
those related to an ensemble of such realizations. 

For a single realization, these quantities are, respec- 
tively: 

(i) The "global" concentration c{t) = N{t)/L^, where 
N{t) is the number of individuals at time t. 

(ii) The distribution of the individual traits P{z, t) and, 
in particular, the time-dependent value of the mean trait. 



^1 ^ Nit) 

(z(t)) = dzzP(zA) = — TT - 
^ " Jo ' N{t) ^ 



(12) 



Another related quantity, currently used by biologists 
[20| . is the mean maladaptation, 



{^iit)) = dz\z 
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(13) 

which represents a measure of the deviation of the pop- 
ulation from the instantaneous optimum imposed by the 
habitat. 

(iii) The distribution of the individual extinction proba- 
bilities n(p, t), that is directly related to P{z, t). Indeed, 
note that Ti{p,t)\dp\ = P{z,t)\dz\, from which one finds, 
using Eqs. ^ and (g]) that: 



n(p,i) 



(1-p) \n'{l-p) 



P{z,t). 



(14) 



For an ensemble of realizations, the main element of 
interest is represented by: 

(iv) The mean extinction time (MET) tcx for a fixed set of 
parameters. Indeed, due to the intrinsic randomness of a 
population's evolution, the extinction time is a stochastic 
variable, that differs from one realization of the dynamics 
to another. As such, one can consider the associated 
probability distribution funtion of the extinction time, 
T(iox) and, in particular, its first moment. 



tc 



(15) 
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A. Constant environment 

Let us consider first a situation that is analogous to 
the one described in Sec. IIV Ai namely a constant envi- 
ronment with if — 0.5, in which evolves a finite popula- 
tion with random initial conditions, which means random 
positions of the individuals on the nodes of the lattice, 
as well as randomly distributed individual traits with 
{z{t = 0)) = 0.5. 

The main question is whether the phase transition 
"extinct-alive" encountered in the MF approach is de- 
stroyed by the fluctuations, or is it still appearing in this 
IBM description. The situation is of course more compli- 
cated, due to the effects of the mutation, i.e., due to the 
existence of the second control parameter represented by 
the amplitude M of the mutation (see Sec. lIIIB]) . A re- 
lated question is, naturally, how is the mutation affecting 
the appearance and location of a phase transition. 
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FIG. 6: IBM approach: The mean extinction time fcx in a 
constant environment as a function of the selection parameter 
S, for different values of the mutation A4. The vertical lines 
indicate the estimated locations of the critical selections for 
each value of A4. The other parameters are c(0) — 0.7, L = 
100, if — 0.5, and the average was taken over 100 realizations. 



1. Evolution of the concentration c(t) 

An indication on a possible phase transition is obvi- 
ously offered by the asymptotic behavior of the concen- 
tration c{t) as a function of the selection parameter S, 
for given initial conditions and a fixed value of A4. In- 
deed, (a) for "large-enough" selection pressures , c{t) is 
decreasing, with random fluctuations, till the absorbing 
c = state; (b) for "small-enough" values of S, in the 
long-time limit the concentration c{t) reaches a regime of 
stationary fluctuations around a non-zero constant mean 
value . These results clearly point to the existence of 
a phase transition "extinct-alive" . Of course, one has to 
be aware that a "true" phase transition can only appear 
in the thermodynamic limit of an infinite population on 
an infinite-size lattice, see also [131 . In particular, due to 
the finiteness of the system, the transition is not abrupt 
(i.e., there is not a clear-cut critical value iSc), and the 
frontier between the basins of attraction of the absorbing 
state c = and of the non-trivial state is also spread 
out. 
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FIG. 7: IBM approach: Snapshots of the probability distri- 
bution function of the trait P{z,t) at two times t = 100 and 
t — 1000. The thin continuous and dashed-dotted lines refer 
to the population without mutation, = 0. The thick line 
corresponds to the population with = 0.1, for which P{z, t) 
reaches after some 10 MCS a stationary profile. The values 
of the other parameters are S = 0.13, c(0) = 0.7, L = 4000 
(for a good statistics), and ip — 0.5. 



2. Behavior of the MET tc^ 

Another good indication about the existence of the 
phase transition "extinct-alive" is offered by the behavior 
of the MET tcx as a function of S. The main results are 
resumed in Fig. [SI that indicates a sudden, very abrupt 
increase in iox with decreasing 5, which is the signature 
of the phase transition. Moreover, Fig. [6] is illustrating 
the way the mutation Ai modifies the location of the 
transition. 



3. Role of the mutation 

It appears from Fig. [5] that increasing mutation ampli- 
tude M is unfavorable to the survival of a population in a 
constant environment and the transition "extinct-alive" 
is shifted to lower values of S. The argument is that a 
larger A4 implies a larger chance for an offspring to be 
ill- adapted to the constant environment, even though the 
parents might be well-adapted. 
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To get a more quantitative picture of the various fac- 
tors that are involved in this result let us consider a given 
population and focus on the behavior of the probability 
distribution function of the trait P{z,t) Let us fix 

the value of the selection - for example, S = 0.165, such 
that, according to Fig. [6l a population without mutation 
(M. = 0) will survive (i.e., its asymptotic regime will cor- 
respond to a constant mean concentration ^ 0). Let 
us compare it with the case of a population with non-zero 
mutation A4 — 0.1, i.e., a population that, again accord- 
ing to Fig. [51 will not survive the SP. Then, according to 
Fig. H for the population without mutation, = 0, the 
distribution of the trait P{z,t) is narrowing towards a 
(5-peak around z — ip — 0.5 by a progressive elimination, 
due to the SP, of the ill-adapted individuals. Moreover, 
the offsprings of well-fitted parents will be also well-fitted, 
since their trait is a pure inheritance equal to the mean 
trait of the parents. On the contrary, for the case of the 
population with Ai = 0.1 one notices that, after a rather 
rapid transient (a few 10s MCS), the distribution of the 
traits remains practically constant. So, although the ill- 
adapted individuals are eliminated in priority by the SP 
(the distribution of the traits tends to be narrowed by 
the selection), at each mating the traits of the offsprings 
are randomized again by the mutation (the distribution 
of the traits tends to be broadened by the mutation). 
The resulting stationary distribution, that has a finite 
width of the order of 7W , is a balance between these two 
tendencies. 

Therefore the population with = 0.1 is at each mo- 
ment more vulnerable than the population with = 0, 
and may thus die even if the population without muta- 
tion is surviving. 

One notices that the evolution of the population, even 
in the case A1 = 0, is much more complex than that 
described by the MF approach. In particular, the dy- 
namics of the concentration and of the distribution of 
the individual traits are strongly inter-related. As long 
as A^ = 0, they are taking place on a comparable time 
scale, but it is difficult to assimilate this time scale with 
the one of the simple relaxational processes that appear 
in the MF description. However, when A4 ^ 0, the evo- 
lution of the concentration and that of P{z, t) take place, 
in general, on different time scales. For example, for the 
situation of a population with A4 = 0.1 considered in 
Fig. [3 , the relaxation of P{z, t) to its asymptotic profile 
takes place during a short transient time of the order of 
10^ MCS, while the relaxation of the concentration c(t) 
to the absorbing state c = is very slow, with a charac- 
teristic time of the order 10^ MCS, that can be estimated 
from the asymptotic part of the plot c{t). 

As can be seen, the mutation is introducing essentially 
new features in the dynamics, even in the simplest con- 
text of a constant environment. In particular, the charac- 
teristics of the population (like concentration, trait and 
individual extinction probabilities) can evolve on differ- 
ent time scales, that are sensitive to both control param- 
eters of selection S and mutation A4. This fundamental 



role of the mutation is not at all accounted for in the MF 
description. 

After discussing these essential qualitative differences 
that appear between the MF and the IBM results, let 
us turn to the problem of the calibration of the param- 
eters corresponding to the two levels of modeling. One 
has to retain that our MF description is not the result 
of a rigorous coarse-graining procedure applied to the 
IBM-model. As such, besides the problem of the corre- 
spondence of the unit-time already mentioned in Sec. [Tll 
the control parameter of selection S that appears in both 
approaches may present some calibration problems, too. 
Moreover the difficulty in establishing a quantitative cor- 
respondence between the two aproaches is actually far 
from being a simple rescaling problem.lt is related to the 
role of the spatial correlations between the individuals 
(like, for example, possible clustering effects), role that 
is completely neglected in the MF approach. Let us con- 
sider just two aspects in order to illustrate this point: (a) 
It seems that these correlations tend to increase the value 
of the asymptotic concentration c'^, however they are 
contributing to reduce the survival chances at increasing 
SP as compared to the MF. (b) The comparison between 
the role that the maximum number of offsprings A'ofi is 
playing in the MF description and in the IBM one indi- 
cates that in both cases the critical value of the selection 
Sc is increasing with A'off . However, the concentration 
of the stationary non-absorbing state is decreasing with 
Aoff in the MF description, while it is increasing with 
A'off in the IBM approach. 



B. Time-periodic changes of the environment 

Let us now turn to the case of a time-dependent habi- 
tat, whose optimum varies in time according to Eq. (jlOp . 

1. Phase diagram: role of selection S and mutation M 

Consider the location of the phase transition point 
"extinct-alive" , i.e., the critical value of the SP for a pop- 
ulation with given initial conditions. This critical point 
depends not only on the amplitude A and period Tp of 
the optimum variation, but also on A4 - an aspect which 
of course cannot be captured by the MF description. One 
is thus led to consider the phase diagram in the plane of 
the characteristic parameters S and A4 , see Fig. [8l 

Several characteristics are emerging from the analysis 
of these critical curves, namely: 

(a) The phase diagram is very little affected by the ac- 
tual value of the optimum (/? of a constant environment. 
Indeed, as illustrated in Fig. [51 the curves for ip = 0.5 
and (f = 0.8 present little differences. 

This result holds provided that the population has 
initially randomly-distributed individual traits. In this 
case, whatever the constant value of the optimum if, 
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FIG. 8: IBM approach: The phase diagram "extinct-alive" in 
the S and M plane, for various periods Tp of the optimum 
oscillation. For comparison we also presented the results for 
a constant environment with <^ = 0.5, respectively </? = 0.8, 
as well as for an abrupt jump of the optimum (at ti = 1000) 
between these two values. The average was done over 100 
runs and the vertical size of the symbols is of the order of the 
estimated error bars. The continuous lines are just a guide 
for the eye. The values of the other parameters are A = 0.3, 
c(0) = 0.7, and L = 100. 



a fraction of the initial population is (more or less) ill- 
adapted and will be progressively eliminated. Also, the 
descendancy will become progressively more and more 
fitted to the environment. As such, in the long run the 
distribution of the trait of the population will become 
peaked around Lp, with a variance of the order of M, 
and the corresponding distribution of the individual ex- 
tinction probabilities will be essentially the same for all 
constant values of ip. The small differences that appear, 
especially for small values oi M., in the critical curves 
for Lp = 0.5 and (p = 0.5 -I- A are due to the way we 
prescribed the trait of a progeny, according to Eq. ([3]). 
When starting from randomly-distributed initial individ- 
ual traits, the heredity part tends to "push" the descen- 
dancy towards a mean trait equal to 0.5; thus the relax- 
ation towards the asymptotic profile of the distribution 
of the individual extinction probabilities will be slower 
for Lp = 0.5 + A than for p = 0.5. This makes the popu- 
lation with (/3 = 0.5 + j4 a little bit more sensitive to the 
SP. 

Also, as already mentioned in the previous section, for 
a constant environment any modification of the status 
quo of the population due to the mutation is "bad" , and 
the critical value of S decreases monotonously with M.. 

(b) Consider now a population with initially equal indi- 
vidual traits Zi = 0.5 and put it in a constant environ- 
ment with ip — 0.5 + A. For small values of the muta- 
tion the randomness in the descendancy traits cannot 



overcome the misfit inherited from the parents. It is thus 
clear that the population will get extinct for smaller val- 
ues of S than in the cases discussed above. This point is 
also reflected in the behavior of the MET. The result for 
Al = is in agreement with the MF case, see Fig. [5] 

However, for large values of M the inherited traits of 
the offspring are greatly changed by mutations and the 
progeny adapts to the constant environment with Lp = 
0.5 + A. The critical curve is approaching the ones for 
the constant optimum discussed at point (a) above. 

For intermediate values of A^, inheritance and muta- 
tion are both important in the offspring traits, and the 
critical curve still lies below the ones corresponding to 
point (a). This is in agreement with the behavior of 
the mean extinction time, as illustrated in Fig. [TO] for 
M = 0.1. 

As such, the population is sensitive to the initial 
distribution of the individual traits only for small values 
of A^, and the memory of the initial state is lost when 
mutation is strong. 

(c) Any variation in the optimum is harmful for the 
system, as indicated by the fact that all the critical 
curves corresponding to a variable optimum lie below 
the one for the constant optimum p ~ 0.5. This is again 
in qualitative agreement with the MF results. 

(dl) Consider the case of a population without mutation, 
A4 — 0. The less favorable situation for such a popu- 
lation seems to be that of an abrupt jump in the opti- 
mum from Lp — 0.5 to p — 0.5 + A, as compared to all 
the cases of smooth oscillation of the optimum around 
Lp = 0.5. This result is in qualitative agreement with the 
one in Ref. [l^, and can be understood as follows: at 
t = ti, when the perturbation of the optimum sets in, for 
a population with A4 — the overhelming majority of 
the individuals (and thus their children) have an individ- 
ual trait Zi = ip = 0.5. An abrupt change in the optimum 
Lp will lead to an abrupt increase in the extinction prob- 
ability of the individuals (per MCS) till its maximum 
possible value p = 1 — exp[— 5/(1 — A)], corresponding 
to the given amplitude A. This leads to a decrease in 
the critical value of the selection S. One can therefore 
conclude that after a sharp change in the living condi- 
tions, such as a large-scale catastrophe, even a moderate 
SP could be lethal for a population lacking the variety 
induced by mutation. 

Moreover, a long period of the optimum variation is 
less favorable to the population survival than a shorter 
one. Indeed, one can argue that in the case of a large 
period the individuals are spending a longer time (in 
MCS) in an unfavorable environment, and thus have 
more chances to die. Recall that the absence of mutation 
leads to the fact that the descendancy is as ill-adapted 
to the environment as the parents are. 

These elements are also illustrated in Fig. [5] for the 
dependence of the mean-extinction time on S. 
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(d2) Consider now a population with an intermediate 
value of the MA, say M « 0.1. One notices that, con- 
trary to the situation (dl) of a population with = 0, 
an increase in the period Tp of the optimum change is 
beneficial for the population and the critical SP is de- 
placed towards larger values. 

Indeed, for the case of an intermediate A4, the heredity 
and the mutation are equally important in establishing 
the trait of the population. As such, if the optimum 
is varying slowly, the population has the time to adapt 
to the instantaneous value of the optimum, and "good 
traits" are transmitted, through heredity, to the descen- 
dancy; moreover, the randomness element due to the mu- 
tation is somehow compensating for the small variation 
of the optimum from one generation to the other. On 
the contrary, if the optimum is varying rapidly from one 
generation to the next one, even if the parents were well- 
adapted at some instantaneous value of tp, they will be- 
come soon ill-adapted and transmit this misfit to their 
progeny through heredity. The mutation element in the 
offspring traits cannot compensate for the too large in- 
herited misfit between the individual trait and the in- 
stantaneous optimum. 

Therefore, a population with an intermediate muta- 
tion amplitude is more vulnerable to rapid variations in 
the optimum. This result is also illustrated in Fig. (TU] 
through the behavior of the mean extinction time of a 
population, as a function of S, for different values of the 
optimum period Tp. 



tion. Through the selection, unfitted individuals tend to 
be eliminated, and the fitted ones survive. If the op- 
timum varies slowly (large Tp), then the "good" traits, 
that will be passed to the next generations, plus a small 
mutation, remain good for some time. A small amount 
of mutation is therefore beneficial. If, however, Tp is too 
short, then quite soon the previously good traits (that a 
too small mutation cannot correct, in order to accomo- 
date them to the evolving value of the optimum !) turn 
out to be bad. A larger effect of the mutation is needed 
to ensure an optimal adaptation of the descendancy to 
the new values of the optimum; however, see (d3), the 
MA should not be "too large" either. Such effects of a 
beneficial range of mutation have already been observed 
in biological systems, see e.g. [32j |. 

This discussion also indicates that an optimal muta- 
tion can only be found for intermediate values of Tp] too 
large or too small rates of changes of the optimum do 
not allow for the above-described effects, as illustrated 
in Fig. [51 The MET can be thus increased by the op- 
timal amplitude of the mutation, for a given period Tp 
of the optimum variation. This result is illustrated in 
Fig- m for the case of Tp = 1000, for which the phase dia- 
gram of Fig. [8] indicates an optimal value of the mutation 
M ~ 0.1. Indeed, for values of A4 above and below this 
opimal value, the MET (for each value of S) is clearly 
lower than the one corresponding to the optimal muta- 
tion. 



(d3) Very large MA is always harmful for the popula- 
tion's survival, whatever the period of change Tp of the 
optimum. Indeed, a large MA changes so much the traits 
that the inheritance of good traits tends to be lost. Each 
generation is therefore composed of too many individuals 
with unfitted traits (since the "cleaning action" that the 
selection operated on the previous generations does not 
lead to a better-fitted offspring). The critical value of S 
is decreasing monotonously with A4 in the whole region 
of high MA. 

(e) The small- Tp curves are, in general, monotonously 
decreasing with increasing M, suggesting that larger 
mutations, combined with rapid changes of the environ- 
ment, are harmful for the population, i.e., reduce its 
resistance to the SP. 

(f) Finally, the most spectacular element of this phase 
diagram is the existence of an optimal value of the mu- 
tation for intermediate values of Tp (at fixed A), i.e., a 
maximum in the critical curve S-A4. For this optimal 
mutation, the critical value of the selection is maximum, 
indicating that the population resists to higher SPs than 
for other values of A^. The value of this optimal mutation 
depends on the period Tp of the optimum variation. 

A crude explanation of this effect would follow these 
lines: Consider first a population with a small MA, that 
brings diversity to the phenotypic pool of the popula- 
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FIG. 9: IBM approach: The mean extintion time tex as a 
function of t iS for different values of A4. The thin vertical 
lines indicate the estimated location of the transition point. 
The average was done over 100 populations with c(0) — 0.7, 
L = 100, A = 0.3, Tp = 1000 and ti = 1000. The continuous 
lines are just a guide for the eye. 



The signature of the phase transition "extinct-alive" 
can be also found in the behavior of the profile of the 
probability distribution function of the extinction time 
T{tox) for an ensemble of populations with fixed initial 
conditions c(0), (z) and given parameters Ai, A, and Tp. 
Indeed, below the critical value of S the distribution of 
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the extinction time is narrow-peaked around tox, while 
above the transition its profile changes, by developing a 
fat tail towards large values of tox- Therefore, above the 
critical point ^cx is no longer corresponding to the most 
representative value of tex- 



2. Role of the period Tp and amplitude A of optimum 
oscillations 

The influence of the period Tp (for a fixed value of A) 
on the critical value of S depends strongly on the value 
of A4. As shown in the previous paragraph, for very 
small or zero values of a decrease in Tp (i.e., slow 
perturbations of the optimum) leads to an increase in 
the critical value of S. 

However, for intermediate values of A4, the role of Tp 
gets reversed, namely increasing Tp (with fixed A) leads 
to an increase in the critical S. This point was discussed 
in the previous paragraph, and it is shown in Fig. [TU] of 
0.1. 



ular, for large MAs Ai the transition point depends only 
weakly on the value of A. 
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FIG. 10: IBM approach: The mean extintion time tex as a 
function of the selection pressure S for M = 0.1 (intermediate 
MA) and different values of Tp. The results are averages done 
over 100 populations with c(0) = 0.7, L = 100, A = 0.3, and 
ti = 1000. 



Finally, for the case of a population with a large MA, 
the MET is less sensitive to the period Tp of the optimum 
oscillation, as well as to the initial distribution of the 
individual traits than in the case of intermediate or small 
values of M . 

The role of the oscillation amplitude A is easy to re- 
sume: increasing A has a destabilizing effect on the sys- 
tem, i.e., it deplaces the critical curve towards smaller 
values of the SP, see Fig. [TTl and it reduces the extinc- 
tion time of the population (all the other parameters, as 
well as the initial conditions being kept the same). This 
effect, however, is more or less pronounced depending on 
the values of the other parameters Tp and A4 . In partic- 
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FIG. 11: IBM approach: Role of the amplitude A of the op- 
timum oscillation on the phase diagram in the (5, A^) plane. 
The average was done over 100 populations with c(0) — 0.7, 
L — 100, ti = 1000, and the period of optimum oscillation is 
Tp — 1000. The continuous lines are just a guide for the eye. 



3. Influence of the system size 

The system size is playing a role on the location of 
the transition point, as well as on its "sharpness", see 
Ref. t27|]. As expected on general backgrounds, the width 
of the transition zone is increasing with decreasing sys- 
tem size. Moreover, small-size populations are, on av- 
erage, more vulnerable than larger-size ones, i.e., it is 
more probable that small-size systems get extinct for 
smaller values of the SP than the analogous large ones. 
This effect is essentially related to the larger amplitude 
of the relative fluctuations in the number of individuals 
in smaller-size systems as compared to larger-size ones; 
indeed, for a given mean concentration, this amplitude 
goes roughly as 1/L. This effect is well-known in biol- 
ogy as the demographic stochasticity effect in population 
extinction, see Refs. [2^[2^ for further comments. 



4. Evolution of the concentration c{t) 

For "large" (i.e., above the transition point) values of 
the SP the concentration c{t) tends to zero, but below 
the transition point the concentration reaches, after a 
transient regime, an oscillatory behavior of period Tp/2, 
as illustrated in Fig. [T^ Of course, due to the inher- 
ent stochastic nature of the dynamics, these oscillations 
are noisy. As in the MF case, the amplitude of these 
oscillations (for given S and M) depends on both A 
and Tp. The oscillations of the concentration are not 
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in phase with the absolute value of the optimum, but 
exhibit a mean lag that depends on the control parame- 
ters A, Tp, S, and Ai of the system. In Fig. [12] we also 
represented the temporal evolution of the mean popula- 
tion trait {z{t)), see Eq. ([T2|). Its oscillations, of period 
Tp, appear only for A4 ^ 0, while for 7W = the mean 
trait settles down to a constant value that is equal the 
mean value of the optimum ((z) = 0.5). Finally, Fig. [T^ 
presents the oscillations (of period Tp/2) of the mean 
maladaptation {n{t)) (Eq. (fT3|) ). and one notices that 
an instantaneous minimum in the concentration corre- 
sponds to a maximum of the maladaptation. Again, the 
characteristics of these curves (amplitude, mean over one 
period, lag as compared to the optimum variation, etc.) 
depend on S and M. 




2000 3000 
t 



5000 



FIG. 12: IBM approach: The oscillations of the concentration 
c{t), mean trait {z{t)} (dashed line), and mean maladaptation 
(/i(i)) for a population that survives in an oscillating habitat. 
The values of the other parameters are M = 0.1, S = 0.13, 
A = 0.3, Tp = 1000, ti = 1000, L = 100, and c(0) = 0.7. For 
comparison we also represented the oscillating optimum 'p{t). 



VI. CONCLUSIONS AND PERSPECTIVES 

In this paper we have considered a model of single- 
species population dynamics on a 2D lattice, and we an- 
alyzed the role of the selection pressure and mutation 
on its behavior in a periodically-changing environment. 
Two levels of description were considered in detail - a 
mean-field one, and an individual-based approach with 
built-in stochasticity. 

The MF approach allowed us to highlight, in a rather 
intuitive manner, the delicate interplay between the dif- 
ferent time-scales of the processes involved in the dy- 
namics of this highly-nonlinear system. In particular, 
the influence of the amplitude and periodicity of the op- 
timum variation on the critical value of the SP (above 



which the population gets extinct) was discussed in de- 
tail. However, this type of MF approach is unsuitable 
for describing another essential element of the popula- 
tion dynamics, namely the stochasticity introduced by 
the mutation in the descendancy traits [39j. 

In order to take into account the effects of the muta- 
tion, one has to appeal to a more refined level of descrip- 
tion like the individual-based one, in which the stochastic 
aspects of the dynamics are fully accounted for. As ex- 
pected, at the IBM level of description the dynamics of 
the population is richer. In particular, the MA is influ- 
encing in a highly nontrivial manner the critical proper- 
ties of the system at the phase transition point between 
"alive" and "extinct". A spectacular effect in the phase 
diagram of the system in the plane of the parameters 5- 
A4 (for fixed values of the period and amplitude of the 
environmental changes) is the existence of an optimal 
MA, i.e., a value of Ai for which the critical SP is max- 
imal. The existence of this optimal MA, that is specific 
to each intermediate value of Tp, is strongly reminiscent 
of the stochastic resonance phenomenon encountered in 
stochastic systems submitted to a periodic deterministic 
perturbation, see e.g. [33|. In a broad sense, this effect 
means that the response of a nonlinear dynamical system 
to a periodic perturbation of one of its parameters can 
be enhanced by an optimal amount of noisy stimulation 
of the system. In our case, the response of the system is 
encoded in the MET, and the noisy ingredient is the ran- 
dom mutation of amplitude M. A more explicit, quan- 
titative mapping of our system on a nonlinear stochastic 
equation for the concentration that exhibits the stochas- 
tic resonance phenomenon is currently under study. 

We also compared the effects of a smooth variation 
of the optimum with those corresponding to an abrupt 
changing of the environment (as, for example, a large- 
scale catastrophic event). We showed that sufficiently 
large mutations can increase dramatically the survival 
chances of the population in case of a catastrophe. 

Variations of the model introduced above, for exam- 
ple using different expressions for the individual fitness, 
or for the dependence of the probability of extinction 
upon the SP, and fitting better some possible experimen- 
tal data, could also be considered. Also, more compli- 
cated than diffusive types of motion, that may pertain 
to a "strategy" of the individuals (e.g., a tendency to 
approach or, on the contrary, to avoid other individuals) 
are currently under study. We believe however that most 
of the qualitative features of our model, resulting from 
the competition of several basic processes, are generic. 
Extensions to ecosystems composed of several compet- 
ing species and submitted to a changing environment is 
presently under investigation. 
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